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Abstract. We propose a quasi-particle model to describe the lattice QCD equation of state for pure SU(3) 
gauge theory in its deconfined state, for T > l.STc. The method involves mapping the interaction part 
of the equation of state to an effective fugacity of otherwise non-interacting quasi-gluons. We find that 
this mapping is exact. Using the quasi-gluon distribution function, we determine the energy density and 
the modified dispersion relation for the single particle energy, in which the trace anomaly is manifest. As 
an application, we first determine the Debye mass, and then the important transport parameters, viz, the 
shear viscosity, 77 and the shear viscosity to entropy density ratio, "q/S. We find that both 77 and rj/S are 
sensitive to the interactions, and that the interactions significantly lower both 77 and 77/5. 
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' 1 Introduction 

, The physics of the non-perturbative domain of QCD, un- 

■ like the perturbative domain, is less understood. The physics 
[ of confinement and quark- hadron transition require a deep 

. understanding of this domain of QCD and is an area of 
' intense research. The best known way to address the non- 
, perturbative QCD is the lattice gauge theory[T]. One of 
' the important goals in lattice QCD is the determination 
of equation of state(EOS) for strongly interacting mat- 
ter. The knowledge of EOS provides a platform to study 
many interesting physical phenomena; in particular, at 
high temperatures, this provides the most realistic EOS 
for the hot and dense matter(QGP) created in heavy ion 
collision experiments. 

An interesting question that arises is whether the lat- 
tice EOS(LEOS) results can be understood in terms of 
quasi-particles which are either free, or at most weakly in- 
teracting. A positive answer to this problem would open 
the doors for developing appropriate effective theories which 
can capture the highly non-trivial results of LEOS with a 
simpler physical picture. In developing such a picture, an 
endeavor of this kind may not be expected to yield satis- 
factory results, if naive parametrizations in terms of quan- 
tities such as the effective mass are employed. Rather, they 
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have to be more in the spirit of the Fermi liquid picture of 
Landau [5] where the energy is a complicated functional 
of the number density. We undertake a similar exercise 
here, for pure gauge theory, and show that such a descrip- 
tion can indeed be obtained in terms of excitations which 
may be looked upon as quasi-gluons - with an effective fu- 
gacity which captures all the interaction effects. We find 
that our agreement with the lattice results is not merely 
qualitative; its deviation is less than one part in a million. 
The method employed here uses and elaborates upon the 
model introduced earlier by us [3ll4l[5] for studying hot 
pQCD EOS. 

As an application of this effective description, we in- 
vestigate LEOS predictions for the viscosity 77, and the 
viscosity to entropy ratio ij/S. These transport parame- 
ters are central to the understanding of the properties of 
QGP which is produced in heavy ion collisions. Indeed, 
recent experimental observations [6] from RHIC strongly 
suggest that QGP created at RHIC behaves like a near 
perfect fluid, having a very low viscosity to entropy ratio, 
77/5 > l/47r miTllHlin]- This impHes that at temperatures 
close to Tc, the quark matter in the QGP phase is strongly 
interacting, and is perhaps in the non-perturbative do- 
main of QCD. These findings arc in accordance with the 
lattice studies which predict that the hot QCD equation of 
state is approximately 10% away from its ideal counter- 
part even at T = 4r. [TnilTIlfni fn] . It should, therefore, 
be natural to employ LEOS to determine the transport 
parameters. However, theoretical studies [MllTS] seek by 
treating the equilibrium state to be that of an ideal gas 
of quarks and gluons. Such an assumption does not seem 
to be justified in view of the lattice results. Consequently, 
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the determination of -q/S requires a revisit where the non- 
ideal nature of the EOS is explicitly incorporated. Fur- 
ther, since its determination is best undertaken in terms 
of a transport equation [H[|15l[TB] . the quasi-gluon picture 
lends itself naturally to undertake that exercise. 

In addition to studying rj/S, we employ the quasi- 
particle picture to extract the Debye mass, via the trans- 
port equation. This exercise allows us to determine the 
value of the phenomenological coupling constant that oc- 
curs in the Yang-Mills and the Vlasov terms in the trans- 
port equation. As an indication of the robustness of the 
model, we arc able to get a complete agreement between 
the lattice and the quasi-particle results. We make a few 
remarks in passing on the implications to heavy quark 
dissociation in QGP. 

Yet another quantity of interest is the bulk viscosity, 
which survives provided that the trace anomaly is non- 
vanishing. We note that since the quasi-particle repre- 
sentation is exact, it automatically reproduces the trace 
anomaly. It is therefore possible to determine, in principle, 
the bulk viscosity as well by using the tarnsport approach. 
It would be of great interest to compare the results so ob- 
tained with those obtained in Refs. [17l[T5] . This study will 
be undertaken separately. 

The paper is organized as follows. In Section 2, we 
introduce the quasi-particle model, and extract the equi- 
librium distribution function from the pure lattice gauge 
theory EOS, and discuss physical meaning of the effective 
fugacity. In Section 3, we study the temperature depen- 
dence of the number density of the quasi-gluons. By plug- 
ging in this expression in the non-abelian Vlasov equa- 
tion, we determine the Debye mass and also the value 
of the phenomenological coupling constant. Using this we 
further estimate the dissociation temperatures for heavy 
quark systems. In Section 4, we determine the temperature 
dependence of gluon quenching parameter, q. We further 
determine the shear viscosity r/ , and the ratio rj/S. We 
do find a very small value for rj/S, as the experiments 
suggest. In fact, we find that it can violate the AdS/CFT 
(KKS) bound ^[S]. In section 5, we present the conclu- 
sions and future prospects. The mathematical details of 
the determination of rj for LEOS has been shown in the 
Appendix. 



2 Quasi-particle model for pure gauge theory 
EOS 

2.1 The effective fugacity 

We now propose a quasi-gluon description of LEOS at high 
temperatures. As mentioned in the introduction, our ap- 
proach is in the spirit of Landau's Fermi liquid theory[2j. 
The quasi-particle description has been introduced by us 
in [2113] . It has been further used in [5] to determine rj and 
r]/S. Yet, the salient features of the model were not fully 
covered in the earlier papers, which we do so here in the 
following. 

The basic idea is to describe the quasi-particles- the 
quasigluons- by a Bose Einstein distribution (see Eq.(IT])). 



As mentioned, the analogy with the ordinary bosons is 
formal. This is so since, as in the Fermi liquid theory, the 
single particle energy of the quasi-gluons (which define the 
distribution) is itself a functional of the number density. 
This functional, which establishes the collective nature of 
the response, is to be determined by employing the lattice 
equation of state. 

We implement the description by writing the distri- 
bution function for otherwise free quasi-gluons in terms 
of an effective fugacity Zg. The effective fugacity contains 
all the interaction effects, and contributes to the energy of 
the gluons in a non-trivial manner. The success of the pre- 
scription is established a posteriori. We obtain an exact 
mapping, and as we show below, the notion of the temper- 
ature dependent effective mass which has been employed 
earlier, is realized only in some limiting situations. We 
caution that the fugacity which we introduce is merely to 
establish the relation between the number density and the 
energy of the quasi-gluons. In short, our problem tanta- 
mounts to determining Zg self consistently from LEOS. 

With the grand canonical distribution function in mind, 
we write the equilibrium distribution for the quasi-gluons 
as 

Zg exp(-/3ep) 



fa = 



(1) 



(l~2,exp(-/3ep)) 

where the quantity {tp = p) would be the energy of the 
gluons in the absence of interactions. The expression for 
Ep^ the energy of the quasi-gluons will be determined be- 
low. It may be noted that Eq([T|) is not the same as the dis- 
tribution function which would follow from a naive adap- 
tation of the Fermi liquid theory. 

On the other hand, the grand canonical partion func- 
tion in terms of the effective fugacity may be written as 
follows, 



HZ) 



V 



(27r)3 



d3pln(l-zgcxp(-/3p)), (2) 



where Vg = 2(iV^ — 1) is the number of degrees of freedom 
for gluons, and V is the volume. As a strategy to deter- 
mine Zg, we Taylor expand the partition function around 
Zg = 1 (ideal gluon gas), and determine the fugacity by 
comparing it with LEOS, order by order. We find that it 
is sufficient to expand upto 0((5z^), where 5z ^ Zg — 1. We 
obtain 

\n{Z)=HZi)+A,^^5z + A2^^{5zf + 0[{5z)% 

(3) 

where ln(Z/) = is the ideal partition function. The 

coefficients Ai and A2 are given in terms of the following 
integrals. 



A2 = 



exp(— u) 



duu 

(1 - exp(-u)) 



2C[3] 



duu' 



exp(-2u) ^l.TT^^ ... 

(l-CXp(-u))2 -3^3 ^™ ^ ' 



It is straight forward to obtain the expression for the 
pressure from the partition function Eq.(I3]) via 



Pgf3V = HZ). 



(5) 
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Fig. 1. (Color online) Behavior of Zg as a function of T/Tc- 



Rewriting the lattice expression for the pressure as Pl = 
P/ + APg (where P/ = is the ideal part of the 

pressure and APg accounts for the non-ideal part of the 
pressure) and match with RHS of Eq.(l5]), we obtain the 
following quadratic equation for 6z, 

A2Sz^ +AiSz-^^^^^ = 0. (6) 

This equation posses two solutions for Sz: 

Of the two roots written above, only the first root is phys- 
ically acceptable. This follows from the requirement that 
\dz\ < 1 and the facts that the discriminant in Eq.Q is 
positive and that the ratio A1/2A2 « —0.92. This choice 
also has the correct limit when APg = {zg = 1). 

We have plotted the effective fugacity (zg) from Eq.i]?]) 
as a function of temperature in Fig.l. From Fig.l, it is easy 
to see that Zg attains its ideal value only asymptotically 
and < Zg < 1. More importantly, it is clear from Fig. 2 
that the quasi-gluon description of LEOS is exact when 
T > 1.5Tc- The deviations arc negligible, being of 0(10^^) 
. This agreement assures the reliability of our results for 
observables such as viscosity with the quasi-gluon picture. 

2.1.1 Physical significance of the effective fugacity 

The physical significance of the effective fugacity intro- 
duced in the present paper is different from that of the 
effective mass employed in [T9 ll20ll2 1( 122'. We show below 
that the effective mass description of jl9..20. 21j emerges 
from our more general framework only as a limiting case. 
In fact, Zg regulates the number density as a function of 
temperature, apart from contributing to the dispersion re- 
lations for the quasi-gluons. We analyze the latter feature 
first. 



Fig. 2. (Color online) Energy density and entropy density as 
a function of temperature determined from the quasi-particle 
model. The corresponding lattice results are also shown. The 
quasi-particle results show almost perfect quantitative agree- 
ment with the lattice results. 



2.1.2 The dispersion relation 

Notwithstanding appearances, the energy of the quasi- 
gluons is not merely given by the relation = p. Rather, 
it should be determined from the fundamental thermody- 
namic relation between the energy density and the parti- 
tion function 

£g = -^df,L7liZg). (8) 

Substituting Eq.(l3]) for the partition function Zg, we ob- 
tain the following interesting expression, 

^a^i^J "^'P^P + ^'^^ H^M- (9) 

The modified dispersion relation for a quasi-gluon reads, 

Ep=p + T^dTlnizg). (10) 

Notably, we sec that the dispersion relation has picked 
up an additional contribution, T^Ot In(zg), which is purely 
temperature dependent. Note that the usual fugacity terms 
for free bosons do not contribute to the dispersion rela- 
tion. The additional term is crucial since it owes its emer- 
gence to the nonvanishing trace anomaly in LEOS. There- 
fore, this additional (purely temperature dependent) scale, 
which gives the non-zero conformal measure, is responsi- 
ble for the bulk viscosity. Interestingly, the presence of 
this scale docs not change the velocity of the gluons, since 
Vg = dpEp. 

We now study the situations under which the effective 
mass prescription would be viable. To that end, we cast 
Eq.([Tni) in the form 

{p + T^dTHzg)r^{p^+m'), (11) 
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which leads to the identification 

m^{T) = 2pT^dTHzg) + {T^dT\n{zg)f. (12) 

The first term in the expression for m?'(T) is linear in 
the momentum apart from being temperature dependent, 
while the second is purely temperature dependent. Thus, 
if Zg is to be realized in terms of an effective mass, the mass 
would have to be momentum dependent. However, in the 
low momentum limit ( ultra soft quasi-gluons), the first 
term becomes subdominant wrt the second term. In this 
particular limit, the effective fugacity can be interpreted 
as a purely temperature dependent effective mass. The 
condition translates to p ^ T^dr In(zg). From this we see 
that RHS of Eq.ini) as T cx). 

Integrating the RHS of Eq.([9]), we obtain the energy 
density as. 

Eg = iPg + Ag, (13) 

whereZig = T^dT\n{zg)J\fg is the trace anomaly and Ng, 
is the quasi-giuon number density, 

A/-. = ^ / d'pf!g. (14) 

We shall study Afg in detail in the next section. 

Before we end this section, we note that effective fugac- 
ity descriptions have been earlier employed in condensed 
matter systems in the last decade. We summarize these 
works briefly. To study the nature of Bose-Einstein (BE) 
condensation transition in interacting Bose gases, a para- 
metric EOS in terms of the effective fugacity has been pro- 
posed by Li et a/|23| , which provides a scheme for explor- 
ing the quantum-statistical nature of the BEG transition 
with interacting gases. Effective fugacity has been used 
for a unitary fermion gas by Chen et a/|24| for studying 
thermodynamics with non-Gaussian correlations. Purely 
as technical tool to distinguish the populations in the con- 
densate state from the others, effective fugacity has been 
employed by Haugerud et aZ,[25] for a BE system of non- 
interacting bosons in a harmonic trap. A similar approach 
has been employed in Refs. pMTTll^ for studying BEG 
with interacting bosons. None of them employs the effec- 
tive dispersion relation which we obtain naturally in this 
work. 



2.1.3 Entropy Density 

The entropy density as a function of temperature can be 
obtained from Eq.®, by employing S = ydr \n{Z). After 
some straightforward manipulation, we get 

S^4^ + ^. (15) 

The first term in the above equation is due to the unmod- 
ified dispersion relation, while the second term is nothing 
but the trace anomaly contribution to the entropy density. 

The behavior of the energy density and that of the 
entropy density are shown in Fig. 2. As expected, they 



match with the lattice results, displaying the viability of 
the quasi-particle model. It will be seen in section 4 that 
the temperature dependence of S will make a substantial 
contribution to the temperature dependence of the ratio, 
?]/S for QGP. 

3 The effective number density and Debye 
mass 

3.1 The number density 

We turn our attention to the number density of the quasi- 
gluons, which need not be the same as that of the inter- 
acting gluons. It is given by 

J^y = -^ J d'pf!gip,zs). (16) 

Using the isotropy of the distribution function and per- 
forming the momentum integral one obtains, 

^f9 = ^VolyLog[3,Zg]. (17) 

Its ideal counter part reads (zg = 1), 

AAx = ^C[3], (18) 

where the function VolyLog[n, Zg] = X^fcLo ■ 
we plot the ratio of the number density of the quasi-gluons 
relative to that of ideal gluons, Rj^f = Ng/Nx, as a func- 
tion of temperature. The ratio is always less than unity 
and approaches the ideal limit asymptotically. 
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Fig. 3. (Color online) Behavior of Rm as a function of T/Tc. 
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3.2 The Debye mass 

The Debye mass is independently determined by lattice 
computations, and as such, there is no need to address it 
again within our model. Yet, we may ask if the knowledge 
of M]j can throw light on the effective coupling constant 
g' which occurs in the transport equation. A determina- 
tion of the effective coupling constant is warranted since 
it contributes to the transport and the thermodynamic 
properties of QGP. In particular, the viscosity - which we 
are interested in this paper - depends on g' (see Eq. (|26| ). 

To that end, wc employ the equilibrium distribution 
obtained in the previous section to and write the permit- 
tivity at zero frequency in the form e(w, k) = 1 + M^/k"^, 
in terms of the Debye mass Mn which is given 

by 

Ml = -2N,{g'f J d^pd,J!^ip, zg), (19) 
which, on an explicit evaluation acquires the form 

ON 

Ml = {g'fp-^^VolyLog[2,Zg]. 



Table 1. Radius for various quarkonia states (in unit of fm) 
taken from Ref. 1361. 



(20) 



Wc now match the Debye mass in Eg. ipO]) with the 
lattice parametrized expression for the Debye mass, A/^ = 
(1.40)g(T)T employed in[3T]- This allows us to identify g' 
to be, 

9' = '-'"'^^^^ . (21) 
y'QVolyLog[2,Zg\ 

Incidentally, the plasma frequency ujp = A/d/VS- 



3.2.1 Dissociation temperatures for quarkonia 

We make a brief digression to estimate the dissociation 
temperature for heavy quarkonia, as predicted by LEOS. 
Recall that a quarkonium state is stable against strong 
decay if the over all mass of the pair of quarks remains 
below the open charm and beauty thresholds. The large 
masses of the charm quark(TOc ^ l.bGev) and the bottom 
quark(mb = A.^GeV) allow a study of their spectroscopy, 
based on the non-rclativistic(NR) potential theorv[5^. One 
favorite choice of the potential in the confined phase is the 
Cornell potential. 



V{r) 



a 
r 



(22) 



where a ~ 0.2GeV^ and a ~ 7r/12 are the phenomcno- 
logical parameters. Employing this form of the potential 
in the NR Schrodinger equation leads to the values of 
the radii(rqq) of various quarkonia states as listed in Ta- 
ble 1 . For the complete list of energy and mass of various 
charmonium and botomonium states, we refer the reader 
to Ref. [36| . Note that these numbers obtained from a NR 
theory give a good account of quarkonium spectroscopy 
(the masses are determined with an less than 1% error for 
all spin averaged states). 

In the QGP phase, due to the the screening of chromo- 
electric field, the quarkonioum bound states survive up to 



qq state 


J/'P 




^' 


r 


Xb 




xi 




rqq in fm 


0.25 


0.36 


0.45 


0.14 


0.22 


0.28 


0.34 


0.33 



Table 2. Dissociation temperature(rd) for various quarkonia 
states (in unit of Tc). Note that Tc is taken to be 0.27 GeV[33]. 
We employ 2-loop expression for the QCD running coupling 
constant at finite temperature|34| 



qq state 




Xc 




r 


Xb 






^11 




1.24 


1.00 


1.00 


2.56 


1.47 


1.07 


1.00 


1.00 



a temperature. One simple way to determine the temper- 
ature at which a particular state dissociates is: whenever 
1/Md < Tqq, where rqq is the rms radius of the state, the 
particular state will not survive in the medium. The equal- 
ity yields the dissociation temperature T^, which wc dis- 
play in table 2, by employing the Cornell potential. These 
estimates are smaller than the other estimates for Td [4jl361 
[571[551[5^HD] ■ and somewhat close to the results obtained 
in Ref. [31[|4T| . But this can perhaps not be taken too se- 
riously since the criterion for determining Td requires re- 
finement. 



4 The shear viscosity 

We now consider the important physical quantity, the 
shear viscosity rj and its ratio to the entropy density, rj/S. 
Determination of rj requires a knowledge of the coUisional 
properties of the medium when it is perturbed away from 
equilibrium. Of the two methods that determine the trans- 
port parameters, viz. the Kubo formula, and the semi- 
classical transport theory, we adopt the latter one in this 
paper, and follow the approach of Asakawa et al. [15j . 

The shear viscosity has two contributions [15], (i) from 
the Vlasov term which captures the long range component 
of the interactions, and (ii) the collision term which mod- 
els the short range component of the interaction. The net 
viscosity is given by l/ij = 1/r/A + l/rjc, where the first 
term gets its contributions from the diffusive Vlasov term 
and the second term gets contribution from the collision 
term. Asakawa, Bass and Miiller [l5j have argued that the 
diffusive Vlasov contributions to the shear viscosity domi- 
nates in the weak coupling limit. We restrict our study to 
determine r]A and tia/S here. We shall drop the subscript 
A hence forth. 

In their work, Asakawa, Bass and Miiller [15] have 
considered the Vlasov term for an ensemble of turbulent 
color fields, but assume that the equilibrium configuration 
is that of an ideal gas of gluons. Such an assumption is 
clearly not admissible while employing LEOS. In an earlier 
paper, we have generalized their work to a perturbatively 
interacting QGP |5]. It was found that the inclusion of 
interactions significantly decreases rj and rj/S. LEOS is 
expected to cause similar significant changes, which we 
estimate now. 
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Fig. 4. (Color online) Gluon quenching parameter g as a 
function of T /Tc. 




Fig. 5. (Color online) The quantity 37r^7?/32(7Vc^ - 1)T^ as 
a function of T /Tc is shown. The upper curve corresponds to 
constant q — IGeV^ / fm and the lower curve is obtained by 
employing the temperature dependence of q. 



The method of obtaining ij has been described at length 
in P^lfTM S], Using the same method, we may write, 



-13 f d^p . 



feq(p), 



(23) 



where A{p) parametrizes the anisotropy in the distribu- 
tion(for details see Ref.[5]). A{p) can be determined by 
the variational procedure from the linearized transport 
equation [I^IITS] with a Vlasov term and a collision term 
computed by Arnold et a/|T4] . It is important to note that 
the work of Asakawa et aZ[l5] is the generalization of the 
work of based on the earlier work of Dupree[43j for the 
non-abelian plasmas. 

The form of A{p) employing ideal EOS has been de- 
termined in [15] in the case of a purely chromo- magnetic 
plasma, and they obtain 



A{p) 



[N^ - \)ElT 



mag : 



(24) 



where g' is the phenomenological coupling and r™°f is the 
magnetic relaxation time. We demonstrate that expres- 
sion for A{p) in the present case is formally the same as 
the one given above. This follows from the fact Zg, which 
captures all the interaction effects in f^q is independent 
of momentum and is purely temperature dependent. Re- 
call that accordingly, the expression for the particle en- 
ergy gets modified to Ep = p + T'^dT[ln{zg)] (see Eq. pO|) ). 
and that the energy density is related to the pressure via 
£ = 3P+A. Keeping these in mind, the same procedure as 
in [15] may be followed which yields Eq. ((24|) . The details 
are given in the appendix. 

The expression for A{p) taking Tm along the light cone [35] 
would then be: 



A{p) = 



(25) 



Here, the lightconc frame is introduced only to relate the 
denominator of the above equation with the gluon quench- 
ing parameter. It must be borne in mind that < E'^ + B'^ > 
is essentially the energy density which must be determined 
by taking only the contributions from the soft modes. 
Thus, 



V = 



157r2C2(g')^ <E^ + B^>T„ 



dp//e,(l + /e,). 

(26) 

Employing the distribution function Eq. ([!]), extracted from 
LEOS, one obtains the following expression for the shear 
viscosity. 



V ■ 



{Nl - If 



l^Tv'^Ncig'f <E^ + B^>T, 



m Jo 



P 



Zg exp(-/3p) 
(1 - Zg exp(-/3p))2 ' 



(27) 

m integral becor 
T^VolyLog[6, Zg]. (28) 



which after performing the momentum integral becomes 

_ W{N^ - 1)2 

~ TT^Ncig'f <E^ + B^>T, 

Note that the shear viscosity for ideal gluons is given by, 

, 16C(6)(7V2 „ y6 



V 



{g'Y <E^ + B'^>T„ 



(29) 



In the above expression, < E^ + B^ > gets contributions 
from the soft modes, and is hence not the standard en- 
ergy density. The Debye mass is a convenient parameter 
to demarcate the soft and the hard modes, whence we per- 
form the momentum integration in Eq.(9) only up to Mjj. 
Denoting the resulting energy density by £^ , we obtain 



1=1 



^ff7[3,1.4g(r)/]), 



1=1 



(30) 



where 7[n, is the lower incomplete gamma function: 
7[n, x] = in- 1)!(1 - exp(-x) Efc=d ^Vfc!)- 
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Since the parameter Tm is unknown, one approach is to 
relate < + B'^ > Tm to the gluon quenching parameter 
g[35] as given by 



q 



3(iV2-l) ' 



(31) 



in writing which we have employed the relation 



£^ , which follows from LEOS. In that case, the expression 
for shear viscosity reads. 



32{N^ - l)T^VolyLog[6, 
^TT^q 



(32) 



The expression for the ratio 77/5 can be obtained by comb- 
ing Eq.(I32D and Eq.fTOl). 

Clearly, what can be determined in this approach un- 
ambiguously is the ratio 



- 1) 



(33) 



Thus, we see that in the approach taken above, the 
problem of determining t] reduces to a determination of 
either the gluon quenching parameter g, or of Tm- There 
are several attempts to determine q. The approach based 
on the twist expansion [44] predicts a value g ~ 1 — 
2Gel/7/m at a temperature Tq ~ (337±10)MeV. The es- 
timation based on the eikonal approximation |45j predicts 
a much larger range of values, between 10— 30GeV'^ / fm. 
Note that the above estimates, which are fitted from the 
data, are by no means precise, and are available only at 
one particular value of temperature. The value of rj also 
inherits the same uncertainty. 

While it is not easy to eliminate the uncertainty in q 
in the above mentioned analyses, we show that it is pos- 
sible to determine its temperature dependence, given its 
value at some temperature, say at Tq. We do so by con- 
sidering the parameter Tm instead. For the plasma under 
consideration, an intuitively appealing way is to relate Tm 
to the plasma frequency o^p, as Tm = Cu!~^, where C is 
the proportionality constant. 

The plasma frequency, cUp for LEOS can be determined 
by employing the quasi-particle model in the expression 
for the chromo-electric susceptibility, in the limit fc ^ 0: 

e{uj,0) = 1 — ^ ■ It is easy to check that ujp — AIjj / \/3. At 
this point we employ the expression for the Debye mass 
determined in the previous section. It is important to note 
that Mo, and hence ujp, are sensitive to the interactions, 
which makes our ansatz plausible. On the other hand, we 
fix C by using the value of q at Tq. Since its value has 
been estimated to be in the range 1 — 2Gev'^ / fm [33], we 
get C « 1.07 — 2.14. This data point completely fixes q 
as a function of temperature. We have not employed the 
other set of values since they are not easy to accommo- 
date within the perturbative frame work which we have 
employed here. We have shown g as a function of temper- 
ature in Fig. 4. It is clear that q has a strong dependence 
on temperature which cannot be ignored in the determi- 
nation of 1]. 




Fig. 6. (Color online) Viscosity to entropy density ratio (^) 
as a function of T/Tc. Note that, we have chosen Tc = .27GeV 

Mi- 



We have plotted the shear viscosity, t] for LEOS as a 
function of temperature in Fig. 5, with the choice C — 
1.07. For comparison we have also shown the values of rj 
when q{T) is assumed to be a constant. The strong de- 
pendence of q on the temperature is clearly reflected in 
the viscosity, with its value getting lowered substantially 
around T = 2.5Tc. We have shown the behavior of rj/S as 




Fig. 7. (Color online) 7?.^ as a function of T /Tc. 



a function of temperature in Fig. 6. It appears that the ra- 
tio may violate the the AdS/CFT boimd, l/47r, marginally 
for T < l.STc. A larger violation of the bound is possible 
at higher temperatures, if one employs eikonal based esti- 
mates for q. 

Let us consider the ratio 7?.^ ~ v/''/ to see how rj for 
LEOS deviates from its ideal counterpart. This ratio is 
model independent to the extent that it does not depend 
on q. In making this statement it is understood that as 
a phenomenological parameter, q is not sensitive to the 
EOS employed [33U3S]- The behavior of TZjj as a function 
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of T /Tc is shown in Fig. 7, from which it is clear that the 
inclusion of interactions significantly decreases the shear 
viscosity. As expected, the ratio 7?.,, asymptotically ap- 
proaches unity. Therefore, the shear viscosity serves as a 
good diagnostic to distinguish the EOS at RHIC. 

To see the extent to which the interactions effect the 
?7/iS, we consider the ratio TZri/s = ^1^1% ■ The behavior of 
Ti-n/s ^ a function of temperature is shown in Fig. 8. From 
Fig. 8 it is easy to see that interactions coming from LEOS 
decrease the ratio rj/S by ~ 35% near l.bTc and « 5% 
near 3Tc. It approaches the corresponding ideal value only 
asymptotically. This crucial observation reinforces the ne- 
cessity of employing realistic equations of state, in par- 
ticular LEOS for determining the transport properties of 
the plasma. Our findings lead to an interesting conclusion 
that both rj and the ratio, rj/ S are good diagnostics as far 
as the effects of interactions are concerned. 




out the effect of trace anomaly and also the collective na- 
ture of these excitations. We have determined the tem- 
perature dependence of the Debye mass which can be ex- 
actly matched with the lattice parametrized Debye mass 
by defining the effective gluon charge in terms of the QCD 
running coupling constant. Employing the quasi-particle 
model, we have determined 77 and the ratio rj/S. In doing 
this, we have determined q as function of temperature for 
LEOS. We have also determined the temperature depen- 
dence of gluon quenching parameter. We find that both 
•q and rj/S for LEOS decrease significantly as compare to 
the ideal EOS. We find that there is a possible violation 
of AdS/CFT bound for rj/S for lattice equation of state. 
It would be of interest to extend this analysis to the full 
QCD EOS and also to study the bulk viscosity. Should 
the quasi-particle model work for full QCD equally well, 
it opens up interesting possibilities of building effective 
theories. 
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6 Appendix 

In this appendix, we show how one determines the anisotropy 
parameter, A{p) for LEOS. We start with the equilibrium 
distribution function f^-q = l/(z~^ exp(/3M.p) — 1), where 
Zg is purely temperature dependent, for the quasi-gluons. 
The action of the drift operator on /eg is given by 



Fig. 8. (Color online) Viscosity to entropy density ratio rel- 
ative to its ideal counterpart as a function of T /Tc- Note that 



this plot is model independent since the ratio, 
dependent of q 



n. 



v/S 



IS mter- 



+/3(«-a)(7/-p)l, (34) 



We further note that we recover the expression for the 
ratio r\IS obtained by Majumder tt al [25] S if we take the 
limit — > 1, as a special case. Our results for rj/S arc at 
a variance with the predictions of [47 p 8 t .51 t ,52 , .5,.46 , .7,.49t 



5 Conclusion and Outlook 

In conclusion, we have developed a quasi-particle model in 
the spirit of Landau's Fermi liquid theory to extract the 
distribution function for gluons from pure gauge theory 
equation of state. We find that the description is exact. 
We show that all the interaction effects can be captured 
in the effective fugacity for gluons. We have determined 
the new dispersion relation for quasi-giuons which brings 



where we recognize that p—df} In(zg) = Ep, is the modified 
dispersion relation. In the local rest frame of the fluid, 
this expression is formally the same as Eq.(6.1) in [15] 
where of course Ep = p. Similarly, the expressions for the 
Debye mass and the continuity equation for the energy 
momentum tensor ( Eqs.(6.3)-(6.7) in [15]) also undergo 
the same modification via the new dispersion relation. 

The final expression for the drift term after imposing 
the energy-momentum conservation is obtained as 



ivd)Ug{p)= feg{l + feg) 



EpT 



(Vw), 



mlE^T^^jEp 
3T^d£/dT 

E,, 



(35) 
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where is the speed of sound. The third term in Eg. ([55)1 
will contribute to the bulk viscosity. To determine the bulk 
viscosity for LEOS, we need to include trace part in the 
ansatz for /i (p, f) [IS] • In this case, the form of the per- 
turbation, r) gets modified as, 



/i 



PtPj 
EpT^ 



{Vu),jAiip) + SijliV ■ u)A2{v) ] . (36) 



The second term in the above expression will generate 
a term proportional to V • it in the force term, and the 
comaprison of this term with the third term in Eq. (|35p 
would lead to the expression for A2 (p) and hence the bulk 
viscosity. Since we are only interested in the shear viscosity 
here, we concentrate on the form of Ai{p) = A(p). 

On the other hand, the force term will have exactly 
the same mathematical form as in [T5](Eq.(6.13)), if we 
consider only the traceless part of velocity gradient in the 
expression for /i(p, r). The same mathematical structure 
of the Force term in this case follows from the isotropy of 
feq{p)- The force term in the case of a purely chromomag- 
netic plasma in the present case will be. 



2 —7nag 



x/eg(l + /egWj(V?/)y (37) 



On comparing of the Force term and the first term in 
the RHS of Eq. ([35|l . we infer that the anisotropy param- 
eter is given by 



AiP) 



{Nl-l)ElT 



(38) 
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